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Abstract 

In reconstructing the top quark momentum from its decay products, one 
must account for extra jets that can result from radiation of gluons. In this pa- 
per we study soft gluon radiation in top production and decay at the Fermilab 
Tevatron. We consider the cases where one or both W bosons decays hadroni- 
cally, i.e., where at least one of the top quarks can be fully reconstructed. We 
show results for gluon distributions at the Tevatron and find that radiation from 
the W decay products contributes substantially in the central region. 



1 Introduction 



The top quark's existence is now well established [|I], Q and top searches have given way 
to top studies. In measuring the top quark's properties — in particular, its mass - 
radiated gluons can play an important role. This is due not only to the relatively high 
probability that ti events are accompanied by additional gluon jets, but also to the fact 
that gluons can be radiated in both the top production and decay processes. This can 
complicate top momentum reconstruction, leading to ambiguities in identifying the 
top quark's decay products, increased systematic uncertainties, or both. These effects 
in turn influence any measurement based on momentum reconstruction, and may also 
bias event selection based on identifying top events via invariant mass cuts. Therefore 
it is necessary to understand the distributions of radiated gluons in top quark events. 

There have been a number of studies of gluon radiation in top quark production 
and decay at e+e" §, |, |, |, 0] and pp colliders g J| |, 0, |11|, 0. In we examined 
additional gluons in the processes qq, gg — > bW + bW~ at the Tevatron pp collider in 
the soft gluon approximation and in [Tl[ we performed the exact calculation. These 
studies included all effects of gluons radiated from the initial state partons, from the 
t and i quarks at the production stage, and from the t and b (and i and b) in the 
t — > bW decay; all effects of the width in the top propagator were also included. 
These studies did not include gluon radiation from hadronically decaying W bosons, 
so that, strictly speaking, they are directly applicable only to dilepton events, in which 
both W's decay leptonically. But in the most useful channel — 'lepton + jets' — only 
one of the W's decays leptonically, and the other decays to quarks. The presence 
of one charged lepton in this channel helps suppress backgrounds without sacrificing 
event rate, while jets from the hadronic W decay allow for direct reconstruction of one 
of the top quarks' decay products. These events can also be accompanied by extra 
gluons. In fact they are more likely to do so than dilepton events, because there are 
more colored particles available to radiate. 

In this paper we extend the study of Ref. || for radiation of soft gluons in top 
production and decay at the Tevatron to include gluon radiation in hadronic W decays, 
W — > qq' . In the next section we present analytic expressions for all contributions to 
the gluon radiation probability in the soft gluon approximation, with a discussion of 
properties of the W decay piece. In section 3, we apply the results of section 2 to 
ti production and decay at the Tevatron and show the resulting gluon distributions. 
Expectations for the LHC are briefly discussed. Section 4 contains our conclusions. 
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2 Soft gluon radiation in tt production and decay: 
general formalism 

At leading order, tt production and decay is given by the subprocess(es) 

a(h) + b(k 2 ) -> t( qi ) + t(q 2 ) -> 6(p0 + + 6(pa) + W~ , (1) 

where afe = qq or for hadron colliders, and the particles' momenta are indicated in 
parentheses. Let us assume that the W + decays hadronically: 

W + ^q(p 3 ) + q'(p,). (2) 

If we are interested in the all jets case, then both Ws decay hadronically and we also 
have 

W~ ^q'(p 5 ) + q(p e ). (3) 

Now consider emission of a gluon with momentum k^. It can come from any of the 
quarks or gluons in the above processes (i.e., any of the particles whose momenta are 
labeled). This includes initial state radiation, gluons emitted by the top (or t) quark 
in either the production or decay process,[] and gluons emitted by the daughter quarks 
of the t and i and the W boson(s). 

Because of the infra-red divergence associated with radiation of gluons, an extra 
jet in a tt event will usually be soft. In the limit of soft gluon emission, several simpli- 
fications occur. First, the gluon momentum fc M does not affect the kinematics, so that 
all of the leading-order kinematic relations still hold.Q Second, the differential cross 
section factorizes into the lowest-order cross section multiplied by a gluon radiation 
probability; following [8], [| we can write 

1 da a s „ _ 

Eg?, (4) 



da dE g d cos 9 g d(f) g An 2 

where dao is the differential cross section for the lowest-order process (i.e., with no 
gluon radiation), E g is the energy of the soft gluon, and ot s is the strong coupling. The 
function T is the sum of 'antenna patterns' of the radiation from the different sources 
listed above. It can be written generically as 

J~ = J~ PROD ~l~ 3~ DEC ~T~ 3~\ INT (5) 



1 Gluons emitted by the top quark in the production and decay processes do not involve separate 
Feynman diagrams; integrating over the t quark's virtuality picks up poles corresponding, respectively, 
to emission in the two stages. See [|J for a discussion. 

2 While this makes the computations simpler, it also means that results of the soft calculation 
cannot be used to study the details of top mass reconstruction, which requires an exact treatment. 
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where J-'prod is the contribution from emission at the ti production stage (including 
initial-state radiation), JF DEC is the contribution from emission off the t and i and their 
decay products at the weak decay stage, and JF INT is the contribution from interferences 
between these emissions. Because we are interested in hadronic W decays, it is useful 
to further decompose ^ DEC explicitly into the contribution from radiation off the tfs 
and 6's, and that from radiation off the quarks from the W decay(s): 

•F DEC = <F DEC.tb 4" J~] DEC,W* (6) 

f DBCiW includes contributions from one or both W decays as appropriate. 

The interference term JF INT contains contributions from the interference (i) between 
radiation in the t decay and that in the t decay and (ii) between radiation in ti 
production and radiation in either decay. (As discussed below, there is no contribution 
to the interference from the W decays because the W boson is a color singlet.) In 
principle JF INT can have substantial effects for gluon energies comparable to the top 
decay width T t 0, However, the observable soft jets that are relevant to the 
pp collider experiments have energies much larger than r t , and so in practice the 
interference terms are numerically small. 

Explicit expressions for T have been presented in |8], ||] for the qq and gg subpro- 
cesses, excluding hadronic W decays. The complete expressions, including in hadronic 
W decays, are 



JFprod = cxkxkz + c 2 k x q\ + c 3 kxq 2 + c A k 2 q x + c 5 k 2 q 2 + c 6 glg 2 + c 7 cf[qi + Csqzcfc , 

^DEc.tb = c 7 [qiq~i + vvp\ - 2gTpi] + c 8 [q^q 2 + P2P2 ~ 2q^p 2 } , 

jf DECiW = Cg p^p 4 + CioP^pe , 

jFint = xi {c2[hpi - faqi] + c 4 [k 2 pi - k 2 qi] + c&[qiPi ~ Oife] + 2c 7 [gipi - 9i?i]} 

+ X2 {c3[hp 2 - hq 2 ] + c 5 [k 2 p 2 - k 2 q 2 ] + c 6 [q~Tp 2 - cf[q 2 ] + 2c&[cfip 2 - ^q 2 ]} 

+ X12 c 6 [prp 2 - qyp 2 - qjpi + qtq 2 ] , (7) 

where the functions pq and Xi are defined by 

p-q 



pq 



p ■ k q ■ k ' 

m 2 r 2 



(qi ■ k) 2 + mfrt 

m 2 t T 2 t (g 1 -kq 2 -k + m 2 t T 2 t ) 

\ l; - \, .. i.V> ...2 1 -.21 \i . l.Y> _2r2j • {k ' 



[( qi -ky + m 2 T 2 } [(q 2 -ky + m 2 T 



In terms of the gluon energy, we have T ~ E~ 2 , and so the cross section (Eq. ([|)) has 
the infra-red behavior da/dE g ~ Eg 1 , as expected. Additional collinear singularities 
arise from the p ■ k denominators when p 2 = 0. 
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The coefficients q depend on the color structure of the hard scattering and can be 
different for the qq and gg processes (C F = 4/3, N = 3): 



coefficient 


gg — > 


gg -> tt 


Cl 


l 

N 


-2C F + 2N + 2Y 


c 2 


2C F - jf 


C F -X-Y 


c 3 


2 
N 


C F + X-Y 


c 4 


2 


C F + X-Y 


c 5 


AT _^ 

2Cf - jf 


C F -X~Y 


c 6 


i 

N 


2Y 


c 7> c 8 


-c F 


-c F 


Cg, Cio 


2C F 


2C F 



Note that the coefficients C7-C10 are involved only in the decay contributions (including 
decay-decay interference) and are therefore independent of the production process. 
The quantities X and Y depend on the gg — > tt subprocess energy and scattering 
angle. Explicit expressions and a discussion can be found in [|J. 

The production, tb decay, and interference contributions to T have been studied in 
Refs. |||[|, to which the reader is referred for further discussion. We shall focus here 
on the W decay contribution. Consider a single hadronically decaying W: W — > gg'. 
Its most important feature, for our purposes, is that the W boson is a color singlet. 
Radiation from its decay products at 0(a s ) does not, therefore, interfere with gluons 
from other sources in the process. Its contribution to the radiation pattern simply 
adds incoherently to the total at this order. 

In fact this contribution is just the familiar quark- ant iquark antenna pattern, whose 
properties are well known and are discussed, e.g., in Ref. |L3|. Its color structure is the 
same as for the decay Z — > qq. Here we remind the reader of a few of its properties. 
Because the q and q from the W decay are taken to be massless, their contribution to T 
depends only on their orientation relative to each other and to the gluon. Furthermore, 
the only other kinematical factor is an overall E~ 2 , which is common to all terms in 
T . Explicitly, 

y"DEC,W = ALyp : — 

2C (1 - cos #34) 

F £2(1-COS03)(1-COS04) ' 

where 6*3, #4, and #34 are, respectively, the angles between the gluon and the quark, 
the gluon and the antiquark, and the quark and the antiquark. Figure 1 shows the 
distribution given in Eq. (||) as a function of gluon angle for several gg angles. We see 
the collinear singularities characteristic of radiation from massless particles. In fact 
most (but not all) of the radiation is concentrated in the immediate vicinity of the two 
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(9) 



quarks, which translates in practice to strong sensitivity to jet separation cuts. We 
also see the string effect — the enhancement of radiation in the region between the 
q and q compared to that outside them. This enhancement can be quite significant 
when the quarks are relatively close together in angle. As the angle increases, the 
enhancement decreases, but does not disappear altogether, and for quarks that are 
back-to-back (Fig. 1(d)) we see the plateau characteristic of initial state radiation, 
where the radiating particles are also back-to-back. 

Finally we note that the preceding discussion of the W decay contribution to the 
gluon radiation is true independent of the source of the W, as a consequence of its being 
a color singlet object. In particular, although we are concerned here with hadronic top 
production at the Fermilab Tevatron collider, hadronic W decays in top production 
at, for example, an e + e~ collider, are taken into account in exactly the same way. At 
the order to which we are working, the W decay antenna simply adds incoherently 
to the contribution from the remaining sources of gluon radiation in any process of 
interest .0 

3 Gluon distributions in ti events at the Tevatron 

3.1 Single hadronic W decay: lepton + jets mode 

The soft gluon radiation pattern described by Eq. (^) can be combined with the lowest- 
order cross section to obtain the full gluon distribution in a particular experimental 
situation. In this section we present the resulting gluon distributions for ti production 
at the Fermilab Tevatron pp collider for the case where one W decays leptonically and 
the other decays hadronically, i.e., the 'lepton + jets' mode. As we have seen in the 
previous section, the contribution from the hadronic W decay adds incoherently to 
the other contributions. The latter correspond to what we would have for the dilepton 
mode, i.e., to the results in Ref. 0. We note that results for the production and 

decay distributions presented below are virtually the same as in where a more 
complete discussion of them can be found. Our emphasis here will be on the new 
features resulting from the W decay contribution. 

We consider the process pp — > ti —>■ blubqq, with gluons generated according to 
Eq. (||), with cio = in Eq. ([?]), for pp center-of-mass energy 1.8 TeV. We take 
m t = 174 GeV/c 2 , M w = 80 GeV/c 2 , and m b = 5 GeV/c 2 and we use MRS (A') parton 
distributions We work at the parton level, with kinematic cuts corresponding 

3 Perturbative interference between radiation off the W and other colored particles in the scattering 
process begins at 0(a 2 s ) with the exchange of two gluons in a color singlet state. This has been 
studied in the context of e + e _ — > W + W~ — > qqqq production jl4|, In general, such interference 
is suppressed by factors of l/N% and (T\y/^) 2 , where u> is a scale characteristic of the gluon energies, 
and is therefore numerically very small ]15[ . 
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roughly to those in the experiments: central production of leptons and quarks, some 
minimum separation of jets from each other and from the lepton, and energy and 
transverse momentum cuts appropriate to detectability of jets but also consistent 
with the soft gluon approximation. The cuts are 





< 


1.5, 


M 


< 


3.5 , 


10 GeV/c < 


Pt 


< 25 GeV/c 




< 


50 GeV , 


10 GeV/c < 


Pt > 




ARig, ARij 


> 


0.5 , 



where rj is the pseudorapidity and (AR) 2 = (A77) 2 + (A0) 2 . In Eq. (pTDj) g represents 
the gluon jet and % and j represent any other detected jet (b, b, or either of the quarks 
from the W decay) or the charged lepton. The r\ g and AR ig cuts eliminate the collinear 
singularities associated with emission from massless particles. 

The resulting distributions in the gluon energy and transverse momentum are 
shown as solid histograms in Fig. 2. The contribution associated with top produc- 
tion is shown as a dotted line, that for JF DEC tb as a dashed line, and that for J-'dec.w as 
a dot-dashed line. (Interference terms are included in the total but are too small to 
appear on the plots.) From the figure we see two features of the W decay contribution 
that will also appear in the rapidity distribution, viz., that it is comparable in both 
size and shape to the remaining decay piece. Recall that the latter comes from the 
decays of both the t and the i, whereas the former involves only the decay of a single 
top. 

The gluon pseudorapidity distribution is shown in Fig. 3; again the W decay piece 
looks very similar to the W piece. We conclude that hadronic W decays contribute 
substantially to gluon radiation in the central rapidity region, more than doubling the 
amount of radiation associated with the top quark decays. 

For purposes of distinguishing the production and decay contributions in the dilep- 
ton case, we also considered in || the behavior of the gluon distribution as a function 
of angular distance ARb g of the gluon from the b quark. In Fig. 4(a) we show the dis- 
tribution in ARtg for very small angles (AR < 0.5). Except for allowing ARb g > 0.01, 
we have retained the cuts specified in Eq. (|I~0"|) . We see that in this region the distribu- 
tion is dominated by the HV decay contribution (dashed histogram) and that as AR 
decreases, the distribution increases and then turns over, displaying the 'dead-cone' 
behavior characteristic of radiation from massive particles ||13|| . In contrast, the distri- 
bution in the angular distance ARwj e t, g between the gluon and one of the (massless) 
quarks from the hadronic W decay, shown in Fig. 4(b), should have a singularity when 
the quark and gluon are collinear. We see that when we reduce the minimum allowed 
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AR W j et g to 0.01, the distribution shows no signs of turning over at small angles. The 
singularity is avoided only by the AR cut. 

While it is true that Fig. 4 shows the difference between emissions from massless 
and massive particles, the effect occurs at angles much too small for the difference to 
be detected experimentally. A less academic view appears in Fig. 5(a) and (b), where 
we look at the same distributions on a larger scale after reimposing the AR > 0.5 
requirement. The distributions look remarkably similar, the only real difference being 
which contribution dominates at small angles: the l ttf contribution (dashed histogram) 
in Fig. 5(a) and the W decay contribution (dot-dashed histogram) in Fig. 5(b). This 
could be used in principle to distinguish between the two decay contributions if one 
knew, from secondary vertex detection, for example, which were the b jets and which 
were those from the W decays (assuming the gluon jet to be the softest one). However 
the results are quite sensitive to both the AR cuts and to fragmentation, which is not 
included in our parton-level analysis. 

A perhaps more promising way to separate out the W decay contribution is to note 
that it is associated with only one of the top decays. If the top quarks are produced 
more or less back-to-back, then some of that separation can be expected to survive 
the decay process. In that case the charged lepton would tend not to be too close to 
the quarks from the other W decay. And since gluons prefer to be near the quarks 
that radiated them, radiation from the hadronic W decay should increase with angular 
distance from the charged lepton. The expected effect is quantified in Fig. 6, where 
we show the gluon distribution as a function of the azimuthal angle difference between 
the lepton and the gluon. We consider the azimuthal angles because the t and t are 
produced back-to-back in the transverse plane. In addition to the cuts of Eq. (|Ttj|), 
we require the angular separation of the b and b to be greater than 45°, and that all 
detected particles (lepton and all jets, including the soft gluon jet) have \rj\ < 1. These 
additional cuts are designed to increase the likelihood that the event is more central and 
the t decay products more widely separated from those of the t. The pseudorapidity 
cut on the gluon has the added advantage that it eliminates a large portion of the 
contribution from production (cf. Fig. 3). In Fig. 6 we see the expected behavior. All 
of the distributions are suppressed below about 30° because of the AR cut. Beyond 
that, the production and HV decay contributions are relatively flat. In contrast, the 
W decay contribution (dot-dashed histogram) starts out at about the same size as the 
other two and increases with increasing angle, reaching a maximum about 50% higher 
than the initial value. All of the increase in the total gluon distribution is due to the 
increase in the W decay contribution. 

Similar effects would be expected in the distribution in azimuthal angle difference 
between the gluon and one of the b quarks, due to the typical angular separation of 
about 7r/2 between the b and W from a single top decay. However, using the lepton is 
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more practical experimentally because there is only one lepton in each event fj and also 
because identifying it and measuring its momentum does not involve the complications 
associated with doing so for jets. 



3.2 Two hadronic W decays: all jets mode 

If both the W + and W~ in a ti event decay hadronically, then there are six jets in the 
final state at leading order. The advantages of this mode are its large branching ratio 
and the fact that, in principle, both the t and t can be fully reconstructed from their 
decay jets. In practice, combinatorics and huge QCD backgrounds make this task 
more difficult. In addition, radiation of a gluon can give rise to an extra jet, further 
complicating matters. 

If the additional gluon is soft, its distribution is given by Eq. (f|) with the coeffi- 
cient cio = 2CV in J~bec,w as given in the table. Now, in addition to the W + decay 
contribution discussed above, we also have one from the hadronic decay of the W~ . 
This new contribution is the same as that from the W + , and it adds incoherently. 
Therefore radiation from the W decays is effectively doubled when both W's decay to 
quarks. This is illustrated in Fig. 7, which shows the gluon rapidity distribution at the 
Tevatron, where both W's contribute to the W decay piece (shown as a dot-dashed 
histogram). The cuts are as in Eq. fllPf ). Their effect here is slightly different from 
above, however, because now there is no charged lepton, and the lepton and neutrino 
(to which no cuts were applied previously) are replaced by two quarks. The jet cuts 
now apply to all of the particles in the final state. This accounts for the fact that 
the W decay contribution in the all jets case is slightly less than twice that for lepton 
+ jets. The transverse momentum, energy, and AR distributions exhibit a similar 
doubling of the W decay contribution, and are omitted for brevity. 



3.3 Expectations for the LHC 

The CERN Large Hadron Collider (LHC) will produce large numbers of top quark 
pairs. Gluon radiation in top production and decay at the LHC differs from that at 
the Tevatron largely because of the difference in collision energies. At the higher- 
energy LHC, the top quarks are more energetic, and the cross section is dominated 
by the gluon-gluon initial state, in contrast to the quark-antiquark annihilation that 
dominates at the Tevatron. The result [[T7|] is a much larger proportion of production- 
stage radiation at the LHC than the Tevatron, with relatively little (on the order 
of 10-20%) radiation occurring in the decays in the dilepton case. We expect the 
effects of radiation from hadronic W decays to be similar to effects at the Tevatron: 

4 Not counting those from semileptonic b decays, which can presumably be eliminated with an 
isolation cut. 
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radiation from a single hadronically decaying W would be comparable to the total 
decay contributions from the tb and tb antennae. Production-stage radiation would 
continue to dominate. We defer a detailed discussion to an exact treatment. 

4 Conclusions 

We have presented results for soft gluon radiation in top quark production and decay 
at the Fermilab Tevatron for the cases relevant for direct reconstruction of the top 
quarks: the lepton + jets and all jets decay modes. Allowing for radiation at 0(a s ) 
from hadronic W decays is straightforward because the If is a color singlet. The 
corresponding contribution simply adds incoherently to the total, with no interference 
with the remaining contributions. For this reason, soft radiation in W decays is 
also easily treated in top production in e + e~ collisions by adding J DBCW to the T 
appropriate to the electron-positron initial state (c/. Eqs. (f§,§D). 

We found that, at the Tevatron, gluons radiated in hadronic W decays contribute 
subtantially to the total amount of radiation in the central region. In particular, gluon 
radiation from each hadronic W decay is comparable both in size and distribution to 
the total tb and tb decay contributions. The W and non-W contributions are thus 
difficult to distinguish. But by considering the distribution in the azimuthal angle 
between the gluon and the charged lepton in the lepton + jets mode we saw that 
some separation is possible. We also note that since the tb and tb decay contribu- 
tions comprise the total decay contribution in the dilepton detection mode, the W 
decay contribution could in principle be isolated by making a direct comparison of jet 
distributions in dilepton and lepton + jets events. 

In attempts to measure the top mass by reconstructing the momenta of its decay 
products, it is very important to take into account the presence of extra jets due to 
gluon radiation. As we have shown in the soft approximation, hadronic W decays 
give rise to significant amounts of gluon radiation. An analysis of the impact of this 
radiation on top momentum reconstruction requires an exact treatment, which we 
plan in future work. 
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Figure Captions 

[1] Single W decay contribution JF DEC (Eq. (§)) to the soft gluon radiation pattern 
in the 6 — <p plane, in arbitrary units, for 6* 34 = (a) 45°, (b) 90°, (c) 135°, and 
(d) 180°. 
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[2] Distributions in (a) the gluon energy and (b) the gluon transverse momentum, 
in tt production, via the subprocesses qq,gg — > bW + bW~, in pp collisions at 
•y/J = 1.8 TeV, with W + — > qq. Contributions from production (dotted lines), 
•7~DEc,tb decay (dashed lines), and .F DECiW (dot-dashed lines) are shown along with 
their totals (solid lines). The cuts are listed in Eq. (|10|). 

[3] Gluon pseudorapidity distributions in tt production, via the subprocesses qq, gg — 
bW + bW~ , in pp collisions at y/s = 1.8 TeV, with W + — > qq. The net distribu- 
tion is shown as a solid line; contributions from production (dotted lines), ^DEctb 
decay (dashed lines), and .F D ec,w (dot-dashed lines) are also shown. The cuts 
are listed in Eq. ([TO). 

[4] Distribution in AR between the gluon and (a) the b quark and (b) the charged 
lepton. Cuts are as in Eq. fllOD , except for allowing (a) ARf, g > 0.01 and (b) 

AR Wj et,g > 0.01. 

[5] Distribution in AR between the gluon and (a) the b quark and (b) the charged 
lepton. The cuts are as in Eq. (|T0|). 

[6] Distribution in the azimuthal angle between the charged lepton and gluon jet. 
The cuts are as in Eq. (TO), except \r)\ < 1 for all detected particles and the 
angle between the b and b must be greater than 45°. 

[7] Gluon pseudorapidity distributions in tt production, via the subprocesses qq, gg — 
bW + bW~ , in pp collisions at \fs = 1.8 TeV, with both W's decaying hadroni- 
cally. The net distribution is shown as a solid line; contributions from production 
(dotted lines), JF DEC tb decay (dashed lines), and f DECjW (dot-dashed lines) are also 
shown. The cuts are listed in Eq. dTOf ). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9605369vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9605369vl 
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This figure "figl-3.png" is available in "png" format from: 
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